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Abstract 

We present a preliminary measurement of the time-dependent CP asymmetry in B° — > K^^Kgtr )^ 
decays based on 431 x 10 6 T(AS) -> decays collected with the BABAR detector at the PEP-II 
asymmetric-energy e + e~ collider at SLAC. In a sample containing 316±22 signal events we measure 
Sk*j = —0.08 ± 0.31 ± 0.05 and C^* 7 = —0.15 ± 0.17 ± 0.03. The uncertainties are statistical and 
systematic, respectively. 
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The radiative decay b — * sj [T] serves as a probe of physics beyond the standard model (SM) [2]. 
In the SM it proceeds at leading order through a loop diagram involving a virtual W~ boson and t 
quark. Because weak interactions involve only left-handed fermions and right-handed anti-fermions, 
the photon in b — > 57 is predominantly left-handed while in b — > 57 it is right-handed. Any 
possible interference between the direct decay B° — » K*° (K^tt )^ and the decay via B° mixing, 
B o ^ B° ^ K*° 7, is suppressed by the small rate of b — >■ s^ R relative to b — » sj L , which is of 
the order m s /mb. SM predictions of the CP asymmetry due to interference between mixing and 
decay are expected to be about —0.02 [3J. As discussed in Ref. [2], left-right symmetric models [I] 
could conceivably produce mixing-induced CP asymmetries larger than 0.5. Some supersymmetric 
models without left-right symmetry [5] have also been shown to permit CP asymmetries of 0(1). 
Because the SM asymmetry is quite small, any significant evidence of a larger asymmetry would 
point to a source beyond the SM. 

Here we report a preliminary updated measurement of the time-dependent CP asymmetry in 
B° -> K* -/ based on 431 x 10 6 T(45) -> BB decays collected by the BABAR detector at the PEP-II 
asymmetric-energy e + e~ collider at SLAC. Previous measurements of this mode were performed 
by BABAR $\ and Belle [7]. 

The BABAR detector is described in detail elsewhere [8j. Most important to this analysis are 
the five-layer, double-sided silicon microstrip detector (SVT), the 40-layer drift chamber, and the 
CsI(Tl) electromagnetic calorimeter (EMC). A detailed Monte Carlo (MC) simulation of signal and 
background processes was performed using the EVTGEN [9 J generator and the GEANT4 package [lOj. 

Time-dependent CP asymmetries are determined using the difference of B° meson proper decay 
times At = t s i g — £ tag , where t s i g is the proper decay time of the signal B (B a i g ) and t tag is that of 
the other B (B tag ). The At for B s i g decaying to a final state / is distributed according to 

e -|At[/r s 

V±{At) = — X [1 ± S f sm(Am d At) =F C f cos(Am d At)] , (1) 

where the upper and lower signs correspond to the tag-side B having flavor B° and B° respectively, 
tb is the B° lifetime, and Arrid is the B° — B° mixing frequency. The Ct coefficient is associated 
with the difference in decay amplitudes for B° — ► / and B° — > /, while the Sf coefficient involves 
interference between the B° — B° mixing and decay amplitudes. We note that direct CP violation 
in B° — » X*°7 decays is predicted to be smaller than 1% in the SM pXJ. The current evidence is 
consistent with this, based on self-tagging B — ► K**y decays |12j . 

We search for B° — > K*°(K ^7r°)7 candidates based on the following criteria, all but one of 
which were used in the previous result. Photon candidates are required to have energy greater than 
30 MeV and must have the expected shower shapes in the EMC. The photon from the B decay, 
also called the primary photon, is required to have an energy between 1.5 and 3.5 GeV in the e + e~ 
center of mass (CM) frame to be consistent with b — ► 57 decays |13j . It must be isolated from 
other charged and neutral clusters in the EMC. Primary photon candidates that form tt° — > 77 
or 77 — > 77 candidates of invariant mass 115 < m 77 < 155 MeV/ c 2 or 470 < m 77 < 620 MeV/ c 2 , 
respectively, when combined with another photon of energy greater than 50 MeV for ir° or 250 MeV 
for 77 are discarded. We select — > 7r + 7r~ candidates from oppositely charged tracks for which the 
confidence level of the vertex fit is greater than 0.1%, the tt + it~ invariant mass is between 487 and 
508 MeV/c 2 (about 3d), and the reconstructed decay length is greater than 5 times its uncertainty. 
We select 7r° — ► 77 candidates with invariant mass between 115 and 155 MeV/c 2 (about 3a) and 
energy greater than 590 MeV in the lab frame. We require the invariant i^g7r° mass m(Kgir°) 
to be within 0.8 — 1.0 GeV/ c 2 , and later use its shape in a maximum likelihood fit. We require 
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I cos#k-*| < 0.9, where 0^* is the angle between the K® and the primary photon in the K^tt rest 
frame. Along with signal candidates, we also reconstruct B + — > K* + (K®ir + )^ candidates subject 
to the same requirements as B candidates, and veto events for which the invariant K®7r + mass is 
within 0.8 - 1.0GeV/c 2 . This is new since the last result, and it removes 12% of the background 
due to non-signal B decays. 

We identify signal decays using two Lorentz- invariant quantities: the energy-substituted mass 
m E s = [\J{s/2 + c 2 p e+e - • TpB) 2 /E 2 e+e _ - \pb\ 2 *) /c 2 and the energy difference AE = E* B - yfs/2, 
where (E e + e - ,cp e + e -) and ps = (Eb,cpb) are the four-momenta of the initial e + e~ system and 
the B candidate, respectively, is the CM energy, and the asterisk denotes the CM frame. We 
require 5.2 < toes < 5.3GeV/c 2 and \AE\ < 250 MeV. To discriminate B decays against continuum 
e + e~ — ► qq (q = u,d,s,c) background we require | cos 0^1 < 0.9, where 9* B is the CM angle 
between the B candidate and the e~ beam direction. We also exploit event topology by requiring 
the ratio of Legendre moments £2/^0 to be less than 0.55, where Lj = J2j \Pj\\ cos^*!*, pj is the 
CM momentum of each particle j not used to reconstruct the B candidate, and 9* is the CM angle 
between p*- and the thrust axis of the reconstructed B candidate. 

After all selection criteria have been applied we find the average candidate multiplicity in events 
with at least one candidate is 1.15. In these cases we select the candidate with tt° mass closest to its 
nominal value [2] , and if there is an ambiguity then we select the one with the K® mass closest to its 
nominal value. We evaluate the selection efficiency using simulated events. We find it is about 16%, 
and combined with the B° -> K*°-f branching fraction, B(B° -> K*°-f) = (4.01 ±0.20) x 10~ 5 [T2] . 
we expect 312 ± 24 signal events. We also expect approximately 35 events originating from non- 
signal B decays (BB background). The rest of the 3677 selected events come from continuum 
background. These two background types are treated separately below. 

For each reconstructed B° — > K*°j candidate we use the remaining tracks in the event to 
determine the decay vertex position and flavor of -Btag- A neural network based on kinematic 
and particle identification information assigns each event to one of seven mutually exclusive tag- 
ging categories |15| . including a category for events in which a flavor tag is not determined. 
The performance of this algorithm is determined using a data sample (-Bfl av sample) of fully- 
reconstructed B° — > Z)W~7r + / p + ja\ decays. The average tagging efficiency is measured to be 
Q = J2c e c (! - 2 w c ) 2 = (31-2 ± 0.3)%, where e c and w° are the efficiencies and mistag probabilities, 
respectively, for events tagged in category c. 

We determine the proper time difference between B s [ g and -Bt ag from the spatial separation 
between their decay vertices in the same way as our previous measurement. The B s [ g vertex is 
reconstructed by combining the trajectory with the knowledge of the average interaction point 
(IP), which is calculated every ten minutes based on two-track events during data-taking. The 
-Stag vertex is reconstructed from the remaining charged tracks in the event [16]. We compute At 
and its uncertainty from a geometric fit [T7j to the T(4S) — ► B°B° system, which takes the IP 
constraint |18| into account. The resolution of At is improved by constraining the average sum of 
the two B decay times (i s j g + t tSLg ) to equal 2t b o, with an uncertainty of \^2t b o. We have verified 
in signal MC that no bias on Sk*^ or Cx*-y results from this procedure. 

The At resolution strongly depends on the number of SVT layers traversed by the pions from 
the Kg. In order for the At information to be useful, we require that each pion have at least 2 hits 
in the SVT, and that cr^t < 2.5 ps and \At\ < 20 ps. About 70% of the events in the data sample 
pass these requirements. The events for which the At information is not used can still contribute 
to the measurement of the Cx*-y parameter as long as they have flavor tagging information. 

We extract signal yields and CP asymmetries using an unbinned maximum likelihood fit to 
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tties; AE, L^/Lq, m(Kg7r°), flavor tag, At, and a At- As stated earlier, we expect a significant con- 
tribution from BB background, so we extract the event yield from this source as well as continuum 
background. The likelihood function is the same one used in the previous version of this analysis, 
and is described in detail in Ref. [18]. We assume that the correlation among the observables 
is small enough that the likelihood function can be constructed as a product of one-dimensional 
probability density functions (PDF). A systematic correction is applied later as a result of this as- 
sumption. All signal PDF parameters are determined using simulated events, except for the flavor 
tag efficiencies, mistag probabilities, and At resolution function parameters, which are determined 
from the i?fl av sample. BB background shapes are also determined from simulation. We use the 
large fraction of background events in the fitted data sample to determine continuum background 
PDF parameters. 

The At PDF for signal events and BB background is obtained from the convolution of Eq. CD 
with a resolution function 1Z(5t = (At — At true ), a At)- The CP asymmetries for the BB background, 
S B g and C B g, are fixed to zero in the fit, and we account for a possible deviation from zero in 
the systematic uncertainty. The resolution function is parameterized as the sum of three gaussian 
distributions [16J. The first two have a nonzero mean proportional to the reconstructed a At, 
accounting for a small bias in At from charm decays of the -Btag- Their width is also proportional 
to a At- The third gaussian is centered at zero with a fixed width of 8 ps. We have verified in 
simulation that the parameters of the resolution function for B° — > K* ^/ events are compatible 
with those obtained from the £>fl av sample. Therefore we use the -Bfl av parameters for better 
precision. We assume that the continuum background contains only prompt decays and find that 
the At distribution is well-described by a resolution function of the same form used by the signal 
PDF. The parameters of the background resolution function are determined in the fit to data. 

Figure Q] shows the background-subtracted distributions for uies and AS for B — > K*°j 
candidates. The background subtraction is performed using the sPlot event weighting technique 
described in Ref. [19j . The curves in the figure represent the signal PDFs used in the fit. Fig- 
ure [2] shows the background-subtracted distributions of At for B°- and B°-tagged events, and the 
asymmetry as a function At. 

We find 316 db 22 signal events with 

S K ^ = -0.08 ± 0.31 ± 0.05 

and 

C K » 7 = -0.15 ±0.17 ±0.03, 

where the first error is statistical and the second systematic. We discuss systematic uncertainties 
below. The statistical uncertainties have been increased beyond what was reported in the fit 
result because we have determined them to be underestimated, using an ensemble of simulated 
experiments in which events were generated from the likelihood PDFs. The scaling factors for 
Sk*^ and Cx* 7 are 1.097 and 1.035 respectively. Because the uncertainty of Ck*-j is larger than 
that obtained from the partial rate asymmetry in self-tagging B — > K*^f decays |12| , we also perform 
the fit with Cx*-y fixed to zero and find 

5k» 7 (C!k» 7 = 0) = -0.07 ± 0.32 ± 0.05. 

The linear correlation coefficient between Sk*^ and Cx* 7 is 0.046. 

We consider several sources of systematic uncertainties related to the level and possible asym- 
metry of the background contribution from other B decays. We evaluate this contribution using 
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Figure 1: Signal and background (inset) distributions for ttt-es (left) and AE (right) obtained with 
the weighting technique described in Ref. [19]. The curves represent the PDFs used in the fit, 
normalized to the fitted yield. 
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Figure 2: Signal distribution for At obtained with the weighting technique described in Ref. [19j, 
with Stag tagged as B° (top) or B° (center), and the asymmetry (bottom). The curves represent 
the PDFs for signal decays in the likelihood fit, normalized to the final fit result. 
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Table 1: B° — > K*°j systematic uncertainties. 



Source 


AS 


AC 


B B Background 


0.029 


0.018 


Bias Uncertainty 


0.034 


0.015 


At Resolution Function 


0.011 


0.018 


Beamspot 


0.004 


0.001 


SVT Alignment 


0.002 


0.001 


PDF Uncertainty 


0.025 


0.010 


DCS Z?tag Decays 


0.001 


0.015 


Total 


0.052 


0.035 



simulated samples of non-signal B — > X s ^f and other B decays. For the former we use the Kagan- 
Neubert model [20] to model the photon energy spectrum and JETSET for the fragmentation of 
the s quark. Since the final state multiplicity predicted by the fragmentation model is significantly 
different from BABAR's measurement |21j . we reweight events according to their multiplicity. From 
these studies we expect to find about 35 BB events, with approximately equal contributions from 
B — > X s ^f decays and other generic B decays. The BB background yield extracted in the fit to 
the data is 22 ± 22 events. We vary S BB and C BB within a conservative range derived from the 
composition of the BB background sample and the CP asymmetry averages reported by the Heavy 
Flavor Averaging Group [12] to assign a systematic uncertainty due to the assumption of zero 
asymmetry in this source. Because the BB yield in data is smaller than expected, we fix it to the 
expected value when we vary its CP asymmetry within the S BB range ±0.41 and the C BB range 
±0.33. We assign uncertainties of 0.029 on Sk*^ and 0.018 on Ck*^ based on these variations. 

In an ensemble of simulated experiments created by generating background events from the 
PDFs and embedding signal events from the full MC simulation, we determined there was bias in 
5^*7 but not in C^* 7 . We apply a correction of ±0.067 to Sk*^, with a systematic uncertainty of 
half the shift. We also assign a systematic uncertainty of 0.015 to Cx*-y based on these tests. 

Systematic effects due to uncertainties in the resolution function, beam spot position, and possi- 
ble SVT misalignment are quantified in the same manner as Ref. [22] . We also include uncertainties 
due to imperfect knowledge of the fixed PDF parameters and shapes used in the fit, amounting to 
0.025 on Sk*^ and 0.010 on C^* 7 . Finally, uncertainties of 0.001 in Sk*-j and 0.015 in Cx* 7 are 
included to account for doubly-Cabibbo-suppressed (DCS) decays of the B tag [23J. The systematic 
uncertainties are summarized in Table [TJ 

In summary we have performed a new preliminary measurement of the time-dependent CP 
asymmetry in B° — > K*°(KgTT )^/ decays. We have found it to be consistent with our previous 
result, as well as with the standard model expectation. 
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